Influence of saving propensity on the power law tail of wealth distribution 
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Some general features of kinetic multi-agent models are reviewed, with particular attention to 
f— I ' the relation between the agent saving propensities and the form of the equilibrium wealth distribu- 

' tion. The effect of a finite cutoff of the saving propensity distribution on the corresponding wealth 

distribution is studied. Various results about kinetic multi-agent models are collected and used to 
construct a realistic wealth distribution with zero limit for small values of wealth, an exponential 
form at intermediate and a power law tail at larger values of wealth. 
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I. INTRODUCTION 



1^3 . 

Multi-agent models of closed economy systems have received considerable attention in recent years due to the 
■ fact that they seem to predict realistic shapes of wealth distribution from very simple underlying d ynamics, basicall y 



pie under ly mg dyna mics, Dasi caiiy 

liaiiiainiiiiiiiiii@ 



Qh' equivalent to kinetic theory of ideal gases in classical statistical mechanics [l|, 

A notable plus-point of these simple models is represented by the ability to reproduce the main features of the empirical 
I ' wealth distributions: a Boltzmann distribution at intermediate values of wealth, and a power law at the highest values 
J> , (see e.g. 0,0,^1, ^31131). The power law form in the tail of the distribution was observed more than a century ago 
00 ' by the economist Vilfredo Pareto |2ij , who found that the wealth of individuals in a stable economy has a cumulative 
C^l \ distribution F{x) oc x~", where a, the Pareto exponent, has usually a value between 1 and 2. 

• In these models N agents interact exchanging a quantity x, which can be interpreted as representing any economic 
entity contributing to the agent wealth, expressed in the same unit of measure, e.g. in monetary units. Depending on 
the parameters of the kinetic model, in particular on the values of the saving propensities {A^} (« = 1, • . • , N) of the 
N agents, the equilibrium wealth distribution can be a simple Boltzmann distribution for Ai = 0ll|,|3,|5j) ^ Gamma 
distribution with a similar exponential tail but a well defined mode Xm > for A^ = Aq > la H , or a distribution 
Q ' with a power law tail for randomly distributed Aj H, . It has been recently recognized |l0lll5j that the observed 
, power law arises from the overlap of Gamma distributions, resulting from (subsets of) agents with similar values of A. 
That is, in systems where saving propensity is distributed according to an arbitrary distribution function /(A), agents 
(~| I relax individually toward Maxwell-Boltzmann distributions, similarly to systems with a global saving propensity Aq, 
<^ i- but with the important difference that in this case the various Gamma distributions with different Aq parameters will 
^ overlap and provide the final (power law) equilibrium distribution. 
• I— I , The aim of the present paper is to further investigate the relation between the saving propensity distribution and 
■ the shape of the final equilibrium wealth distribution, with particular attention to reproduce a realistic distribution. 
^ [ In Sec. ^we recall the main features of kinetic multi-agent models, while in Sec. Illll we consider how the equilibrium 
- - ' distribution is affected by a particular choice of the parameters of the saving propensity distribution and provide some 
examples. Results are summarized in Sec. IIVI 
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FIG. 1: Linear (left) and semi-log (right) plots of the probability density of wealth x from numerical simulations (dots) for 
various values of the global saving propensity A in the interval (0,1), compared with the theoretical curves (continuous lines) 
defined by Eqs. JHJ, and IIUI . The curve for A = is the Boltzmann distribution. 



II. KINETIC MULTI-AGENT MODELS 



In kinetic multi-agent models N agents interact with each other through a pair interaction ~ but this is only one 
of the many possibilities - exchanging a quantity x, generally referred to as "wealth" in the following. Agents are 
characterized by their current wealths {x^}, i = 1,2, ... ,7V and, possibly, by some parameters, such as the saving 
propensity A^. The evolution of the system is then carried out in the following way. At every time step two agents i 
and j are extracted randomly and an amount of wealth Ax is exchanged between them, 

X ^ — X X ^ 

x'j = + Ax . (1) 

It can be noticed that in this way the quantity x is conserved during the single transactions, x[ + x'j = Xi + Xj, where 
x[ and x'j are the agent wealths after the transaction has taken place. 



1. The basic model 



In the basic versions of the model the quantity x represents money and Ax the money exchanged, assumed to have 
a constant value 0, S 01 j 



Ax = Axo , (2) 



or to be proportional to the initial values 0, 



Ax = exi — eXj , (3) 

where e is a random number uniformly distributed between and 1 and e = 1 — e. The form of Ax in Eq. Q represents 
a random reshuffling of the wealths of the two agents ^ , since Eq. can in this case be rewritten as 

X^ ^{xi ^~ Xj ) , 

x'j = e{xi +Xj) . (4) 

These dynamics rules, together with the constraint that transactions can take place only if x^ > and x^ > 0, lead 
to an equilibrium state characterized by an exponential Boltzmann distribution, 

/(x) = (x)"^exp(-x/ (x)) , (5) 

where the effective temperature T\ of the system is just the average wealth (see curve A = in Fig. Despite 
its intrinsic simplicity, the basic model has the merit of having shown that economic interactions can be modeled in 
terms of simple statistical mechanisms leading to corresponding universal statistical laws. 
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TABLE I: Analogy between kinetic and multi-agent model 
Kinetic model Economy model 

variable A'^kinetic energy 2:=wealtli 

units A'' particles A'^ agents 

interaction collisions trades 

dimension integer D real number D\ 

equipartition theorem fceT = 2 {K) /D Ta = 2 {x) /Da 

reduced variable ^ = K/ksT ^ = x/Tx 

distribution /(^ = 70/2(6 /(C) = 1d^/2{0 



2. Models with a global saving propensity 

A first generalization toward a more realistic model is based on the introduction of a saving criterion. Agents save 
a fraction A (the saving propensity, with < A < 1) before entering a trade and only exchange the remaining fraction 
(1 — A) of their wealth, 



— Xxi + e(l — X){xi + Xj) , 

= Xxj + e{l - \){xi + Xj) , (6) 



corresponding to a Aa; in Eq. (QJ given by 



Ax ^ {1 - X)[exi - exj] . (7) 
The corresponding equilibrium distribution is well fitted by the gamma distribution [l^ Il3j| 

fiO = p(^j2y e^^'-' exp(-0 ^ 7Da/2(0 , (8) 
as shown in Fig. ^ Here the dimensionless variable 

is just the variable x rescaled with respect to the effective temperature T\ and 

D\ 3A _ 1 + 2A 

^i^W. (10) 

The parameter D\ plays the role of an effective dimension, since the Gamma distribution 7n(C) given by Eq. (|SJ) is 
identical to the Maxwell-Boltzmann distribution of kinetic energy for a system of molecules at temperature Tx in Dx 
dimensions (of course only for integer or half- integer values of n = Da/2) In further support of this analogy, 

it is worth noting that Tx and Dx are related to each other through an "equipartition theorem" , 

, , D\Tx , . 

= • (11) 

The equivalence between kinetic theory and closed economy models , su ggested by the basic version of the kinetic 
multi-agent models Q, |^ |^ |^ , can thus be extended to values A > , as summarized in Table ^ 

While A varies between and 1, the effective dimension Dx increases monotonically between 2 and 00. In fact in a 
higher number of dimensions the fraction of kinetic energy exchanged between particles during a collision is smaller. 
At the same time, the market temperature Tx decreases with increasing A, signaling smaller fluctuations of x during 
trades, consistently with the presence of a saving criterion, i.e. A > 0. One can notice that = (1 — A) (x) + 2A) w 
(1 — A) (x) is on average the amount of wealth exchanged during an interaction between agents, see Eqs. (jB)). 
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FIG. 2: Log-log plot of wealth distribution (dotted line) resolved into partial distributions (continuous curves) corresponding 
to agents with a saving propensity belonging to a sub-interval of the A range (0,1). Left: Partial distributions from the 10 
intervals of width AA = 0.1 of the A range (0,1). Center: Partial distributions from the interval A = (0.9, 1.0) further resolved 
into partial distributions from sub-intervals of width AA = 0.01. Right: Partial distribution from the interval A — (0.99, 1.00) 
- not shown in the central figure - resolved into partial distributions from intervals of width AA — 0.001. Peaks are due to the 
power law breaking down as the distance between consecutive partial distributions becomes comparable with their width. 



3. Models with a continuous distributions of saving propensity 

As a further generalization, various investigations concerned models in which agents have realistically been diversi- 
fied from each other by assigning them different saving propensities A^ |^ Hol jlll 111 . ITsl . In particular, uniformly 
distributed Xi in the interval (0,1) have been studied numerically in Refs. |d. [ici| . This model is described by the 
trading rule 

x[ = X^Xi + e[(l - Xi)x^ + (1 - , 

x'j = Xxj +e[{l~ X^)xi + {1- Xj)xj], (12) 

or, equivalently, by a Ax - as defined in Eq. - given by 

Ax = e(l — Xi)xi — e(l — Xj)xj . (13) 

One of the main features of this model, which is supported by theoretical considerations 0,0,112, is that the wealth 
distribution exhibits a robust power law at large values of x, 

fix)^x-^-\ (14) 

with a Pareto exponent a = 1 largely independent of the details of the Ai-distribution. As remarked in Ref. T^, the 
wealth distribution of the single agents are not of a power law type but have a well defined mode and an exponential 
tail, similarly to the case a global saving propensity Aq. The power law actually results from the overlap of these partial 
distributions corresponding to the various A's, which are Gamma distributions, whose average value is proportional 
to 1/(1 — A) and thus extend to very large values of x J^]. These results are also in agreement with theoretical 
approaches to kinetic multi-agent models [Uj llJi ■ This phenomenon is illustrated in Fig.|5| 

Arbitrarily small (random) irregularities in the distance between two consecutive values of A close enough to 1, in a 
uniform distribution of saving propensity, are amplified in the wealth distribution as a consequence of the correlation 
between average wealth and saving propensity, resulting in isolated peaks in the wealth distribution |15| . This is 
shown by the simple example in Fig.|21 where two distributions, in principle equivalent to each other since associated 
to a uniform A-distribution in (0,1), look actually quite different: The first distribution (left) has been obtained by 
randomly extracting the values of A with a random number generator in the interval (0,1). It can be noticed that 
the corresponding equilibrium wealth distribution is more irregular than that obtained from the second distribution 
(right), in which the values Xi were set to be equidistant from each other in the interval (0,1) by defining them as 
Xi = i/N. The deterministic and random uniform distributions are equivalent to each other in principle but not 
in practice (within numerical simulations), where a finite number of agents is necessarily used. The reason is that 
the A values extracted randomly present fluctuations and therefore wider intervals between neighbor values which 
are amplified in the final wealth distribution. Since in these numerical simulations one tries to mimic continuous 
distributions by use of the smallest possible number of variables, it is convenient to avoid the irregular fluctuations 
present in randomly extracted sets of numbers and use a deterministically extracted sets {A.^} of saving propensities. 
This can be achieved easily with the method prescribed in Appendix A. 
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FIG. 3; Reproduced from Ref. ITsIl : Wealth distribution of a system of A'^ = 10'' agents after 10^^ trades for two uniform 
distributions of A: A randomly generated A-distribution (left) produces a wealth distribution which is more irregular than a 
deterministic uniform A distribution (right) in which A; — i/N. 



III. CONSTRUCTION OF A REALISTIC MODEL 



Here we study how some aspects of the A-distribution /(A) influence the equilibrium form of the wealth distribution 
/(a;), i.e. its shape at smaller and the tail at larger values of and in particular under which conditions an exponential 
and a power law can appear in different ranges of the same distribution. 



A. Wealth distribution at small and intermediate values of wealth 



The equilibrium distribution of the basic model is the simple exponential function in Eq. |(SJ). Such a form of 
distribution decreases monotonously with x and does not have rich agents nor a power law tail, a point dealt with in 
greater detail in the next section. In the small x limit, the exponential distribution is f{x ^ 0) > which implies 
that many agents have a wealth x « 0. In fact the mode of the distribution is Xm = and the fraction of agents 
outside a given interval (0, x) - which is just the upper cumulative distribution function - has a pure exponential 
form, F{x) = cxp(— a;/ {x)). Real data about wealth and income distributions, on the other hand, show that wealth 
distribution functions have a mode Xm > The introduction of a (global) saving propensity 

A > solves this problem since it leads to an equilibrium Gamma distribution 0, : which has a mode Xm > 
and a zero limit for a; — > 0, see Fig.^ This functional form has been shown to interpolate well real data about income 
distributions at small and intermediate values 



B. The tail of the wealth distribution 



The tail of wealth distributions is known to follow a power law with Pareto exponent between 1 and 2, depending on 
the sample analyzed. The model under consideration, when saving propensities are continuously distributed, predicts 
a power law tail in f{x), despite with a lower Pareto exponent a = 1, a feature which has been shown to be very 
robust and independent of the details of /(A). 

Both numerical and theoretical analyzes of kinetic multi-agent models show that agents with large values of A's 
(i.e. A close to 1) give a major contribution to the power law tail. This is illustrated e.g. by the fact that when 
the power law is decomposed into partial distributions of agents within a given interval of the saving propensity, the 
partial distribution corresponding to the interval with the highest A is in turn a power law, while the distributions 
corresponding to lower values of A are localized and have an exponential tail, as shown in Fig. [5] However, they sum 
up to give a power law at lower values of x. All this suggests that the crucial factor for having a power law extending 
beyond a certain value x is the highest A present in the sample, that is the cutoff of the A-distribution. 

Thus, rather than varying the functional form of /(A), the influence of the cut-off Am of the A-distribution ~ which is 
a parameter characterizing numerical simulations as well as real systems - has been analyzed. A uniform deterministic 
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FIG. 4: Wealth distribution obtained from a uniform saving propensities distributions /(A) defined by Eq. I)f 5^ in a system of 
10"" agents for various values of the cutoff Am. Curves from left to right correspond to increasing Am. 



distribution of saving propensity for the N agents in the interval (0, Am), has been generated through the formula 

A. = (^-^^Am, i = l,...,iV, Am<1, (15) 

as described in greater detail in Appendix A. 

In fact we found that varying the cutoff Am influences in turn the cut-off of the wealth distribution f{x) and the 
shape of the distribution at small x - but not the shape of the tail which remains a power law with exponent a = 1. 
Decreasing Am has the effect to decrease the interval of wealth x in which the power law appears, until it eventually 
disappears for Am ~ 0.92. Results are shown in Fig. ^ where the various curves represent the distribution functions 
obtained for some values of the cutoff Am chosen in the interval Am = (0.9, 0.9999) in a system of 10^ agents. Curves 
from left to right correspond to increasing values of cutoff. The transition from an exponential to a power law form 
of the wealth distribution, as the cut-off Am decreases, takes place by a shrinking of the power law interval, rather 
than as a change of the functional form of the tail. 

As a final remark it is to be noted that the cutoff of the A-distribution is naturally linked to that of the x-distribution, 
as a consequence of the correlation existing between average wealth {x) and saving propensity 15] in this model, 

{xi) (1 — Xi) ~ const . (16) 

Here the constant on the right hand side of the equation is the same number for all the agents in the system. This 
relation clearly shows that the highest average wealth is determined in turn by the highest A^. 



C. Superposing an exponential form at intermediate values and a power law tail 

In real wealth distributions, an exponential form at intermediate values of wealth is known to coexist with a power 
law tail at larger values jl^- The power law is mainly due to a small percentage of population, of the order of a few 
per cent, while the majority of the population with smaller average wealth give rise to the exponential part of the 
distribution. In this section we try to construct a realistic example of such a type of wealth distribution by collecting 
some of the results obtained so far: 

• A global saving propensity Ao is associated to an equilibrium Gamma distribution, which always has an expo- 
nential tail. 

• A set of agents with a continuous A-distribution produces a power law in the equilibrium wealth distribution. 
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FIG. 5: Numerical equilibrium wealth distribution of a population of 10^ agents. One per cent (1000) of the agents have 
uniformly distributed saving propensities in the interval (0,1), while the rest have A = 0. Left (semi-log scale): Small x part 
of the numerical distribution (dots) compared with an exponential function oc exp(— 3a;/5)(dotted line). Right (log scale): 
Numerical distribution (continuous line) compared with a power law ex (dashed line) and the same exponential function of 
the left figure (dotted line). Peaks at high x are due to the finite resolution in A, see Fig. |21for an explanation of their origin. 



• The average wealth (xi) of an agent and the corresponding saving propensity are linked to each other through 
Eq. H16() , which imphes that agents with high A « 1 contribute to the large- a; part of the distribution. 

It is then natural to ask if a suitable A-distribution may lead to the desired equilibrium wealth distribution. To answer 
this question we have constructed a hybrid A-distribution - on the base of the results listed above and a very similar 
prescription mentioned in Ref. 'sj- in the following way: 

• A small fraction of agents po with saving propensities A^ uniformly distributed in the interval (0,1) according to 
Eq. KHt . 

• The remaining fraction 1 — po with a constant value of the saving propensity Aq. 

The corresponding distribution for pq — 0.01 (1 per cent) and Aq = is shown in Fig. [Sj both in the small x-scale, 
where the distribution has an exponential shape, and in the long x-range, where the power law with exponent —2, 
which characterizes this type of multi-agent model, is observed. 

It is noteworthy that the coexistence of an exponential and a power law tail is possible only for small values oi po, 
in agreement with the fact that it is a small percentage of the population in real systems that is responsible of the 
power law form of wealth distribution at large values of wealths. For larger values of pq the exponential part shrinks 
and the power law dominates. This effect is in a sense contrary to that considered in Sec. Ill 3l where decreasing the 
cutoff of the A-distribution induced a shrinking of the power law range. 

It may be noticed that, due to the choice Aq = for that part of agents with a constant saving propensity, the 
distribution in Fig. Estill has a mode Xm = 0. However, one recovers a distribution with a well defined mode x^ > 
as soon as one chooses a Aq ^ 0. The distribution in Fig. corresponds to a Ao = 0.2 for 99% of the agents and a 
uniform A-distribution for the remaining agents. 

D. Meaning of the saving propensity 

The central role of the saving propensity A - or risk aversion as referred to in Ref. Q - for the considerations made 
above is evident. However, it is to be remarked that the relation between saving propensity Xi of an agent and the 
corresponding average wealth (xi) should not considered to be of a cause-effect type. It is true that in the present 
model the A^'s are fixed parameters, so that the natural dynamical interpretation is that the saving propensity A 
determines the final average wealth. However, in a real situation the value of A itself may vary according to various 
factors, e.g. the wealth itself: a high average wealth probably puts the agent in a situation which allows to carry on 
trades investing the same amount of wealth while saving more respect to agents with smaller wealths. Therefore, the 
model contains in its very dynamics a positive correlation between A and {x) supported by real data |2^ but leaves 
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FIG. 6: Numerical equilibrium wealth distribution of a population of 10^ agents, in which 10* agents (1%) have uniformly 
distributed saving propensities in the interval (0,1), while the remaining 9.9 X lO'^ agents (99%) have A = 0.2. Respect to Fig.|Sl 
the distribution has a mode Xn^ > 0. Left (semi-log scale): Small scale part of the numerical distribution (dots-continuous line) 
compared with a function oc exp(— 3x/2) (dotted line). Right (log scale): Numerical distribution (continuous line) compared 
with a power law oc (dashed line) and with the same exponential function of the left figure (dotted line). 



the question of the actual dynamical relation between them to a more detailed microscopic analysis. Multi-agents 
models like that considered here describe flux of wealth on a mesoscopic level, i.e. on a coarse grained scale in time 
or wealth, rather than reflecting the single agent strategy to save or reduce risks. 

IV. CONCLUSIONS 

We have shown that within the framework of kinetic multi-agent models it is possible to obtain realistic wealth 
distributions f{x) characterized by a zero limit for small x, and the coexistence of an exponential form at intermediate 
and power law tail at larger values of x. In agreement with observations on real systems, this is possible only if the 
percentage of rich agents does not exceed a critical threshold of the order of 1 per cent. Also, the model naturally 
produces a positive correlation between average wealth (x) and saving propensity A exhibited in real data samples. 
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APPENDIX A: EXTRACTION OF A VARIABLE A WITH CUMULATIVE DISTRIBUTION F(A) 

It is possible to define a sequence of N numbers Ai, i = 1,...,A'', which becomes distributed according to an 
arbitrary distribution function /(A) = dF(X)/d\ in the continuous limit {N — > cx)), in at least two ways, randomly 
or deterministically. The two methods are equivalent to each other only in the continuous limit, while in numerical 
simulations a finite N is necessarily employed and they may provide different results. As discussed in Sec. Ill 31 in 
some cases it may be preferable to have a regular, rather than a randomly extracted sequence. 

• Random extraction. A generator of random numbers 4>, < ((> < 1, uniformly distributed between and 1, 
can be employed to extract a set of numbers Ai distributed in the continuous limit according to an arbitrary 
cumulative distribution function F{X), with F{0) ~ and -F(l) = 1. The cumulative distribution function for 
the random variable (j) is simply F{(j)) = (p and dF = d(j) is the (constant) probability to extract the next random 
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number between (/) and (j) + dcf). The algorithm is based on the identity d(j) = dF = f{X)dX, which shows that if 
values Fi are extracted randomly and uniformly in the interval (0,1), then the corresponding values obtained 
by inverting F — F{Xi) will be distributed with probability density /(A). 

• Deterministic extraction The same result can be obtained by a deterministic assignment of the values A^ which 
does not make use of random number generators. If the sequence {Xi} is assumed to be labeled in increasing 
order, i.e. 0<Ai<A2<---<AAr<l, then the function of i 

AW = A,, (Al) 

increases monotonously with i and it is possible to invert it to express j as a function of A^ to define the function 

^(A.) = ^, (A2) 

which represents the fraction of agents with saving propensity less or equal to A^: F{X) is just the (lower) 
cumulative distribution function and as such < F{X) < 1 for every A, F{X ^ 0) — > 0, and F{X 1) ^ 1. 
For instance the cumulative distribution function of a uniformly distributed variable A in the interval A £ (0, 1) 
is just the linear function F{X) = A, with < A < 1. Then Eq. (|A2|I provides the corresponding deterministic 
sequence as 

A. z = l,...,7V. (A3) 

If there is an upper cutoff Am in the distribution, the equation is modified as in H15|) . In the general case of a 
given cumulative function F{X), it is sufficient to invert Eq. I|A2|1 to obtain the sequence in the form A^ = X{i/N), 
i = 1, . . . , N, where A(. . . ) is the inverse function of F{. . . ). The values Xi thus obtained will be distributed in 
the continuous limit with a probability distribution /(A) — dF{X)/dX. 
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